Femtosecond formation of collective modes due to meanfield fluctuations 
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Starting from a quantum kinetic equation including the mean field and a conserving relaxation- 
time approximation we derive an analytic formula which describes the time dependence of the 
dielectric function in a plasma created by a short intense laser pulse. This formula reproduces uni- 
versal features of the formation of collective modes seen in recent experimental data of femtosecond 
spectroscopy. The presented formula offers a tremendous simplification for the description of the 
formation of quasiparticle features in interacting systems. Numerical demanding treatments can 
now be focused on effects beyond these gross features found here to be describable analytically. 
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The last ten years have been characterized by an enor- 
mous activity about ultrafast excitations in semiconduc- 
tors, clusters, or plasmas by ultrashort laser pulses. The 
femtosecond spectroscopy has opened the exciting pos- 
sibility to observe directly the formation of collective 
modes and quasiparticles in interacting many-body sys- 
tems. This formation is reflected in the time dependence 
of the dielectric function 0, 0, Q or terahertz emission 
For an overview over theoretical and experimental 
work see 0, @. Such ultrafast excitations in semicon- 
ductors have been satisfactorily described by calculating 
nonequilibrium Green's functions This approach 

allows one to describe the formation of collective modes 
0, an d even exciton population inversions [1(1 ] . 

The experimental data in semiconductors like GaAs 
or InP 0] reveal similar features. These features are 
explained by several numerically demanding calculations. 
Since some common features are robust, i.e., independent 
of the actual used material parameters, it should be pos- 
sible to describe them by a simple theory. The origin of 
such robust features likely rests in the short- or transient- 
time behavior itself. At short times higher-order correla- 
tions have no time yet to develop, therefore the dynamics 
is controlled exclusively by mean- field forces. 

Based on the mean-field character of the short-time 
evolution we intend to derive a time-dependent response 
function from the mean-field linear response. We will 
start from a conserving relaxation-time approximation 
which dates back to an idea of Mermin (lHll^. Our aim 
is a simple analytic formula suitable for fits of experi- 
mental data. 

The electron motion is controlled by the kinetic energy 
£ , the external perturbation l/ ext and the induced mean- 
field potential V ind . The corresponding kinetic equation 
for the one-particle reduced density matrix p reads 
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with the relaxation time r and a local equilibrium density 



matrix p determined in the following. We will calculate 
the response in the momentum representation 



f(p,q,t) = (p+ - 
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where \p) is an eigenstate of momentum p. For interpre- 
tations it is more convenient to transform the reduced 
density matrix to the Wigner distribution in phase space 



(3) 



where R is the spatial coordinate. 

The relaxation process in Eq. (JIJ tends to establish 
a local equilibrium. Following Mermin, the local equi- 
librium can be characterized by a local variation of the 
chemical potential. Up to linear order it can be written 
as 



P ' 



MQ,*) (4) 



with the Fermi-Dirac distribution f p — /fd (£p) and the 
kinetic energy |) = e p S(q). The delta 

function S(q) expresses that the corresponding terms are 
diagonal in momentum or translational invariant. The 
local deviation Sp from the global chemical potential is 
selected so that the density of particles is conserved at 
each point and time instant. Momentum and energy con- 
servations can be included, too, leading to slightly more 
complicated formulas |13l 1 1 - 11 ] . 

Let us now specify Sp from the condition of density 
conservation which requires that the local equilibrium 
distribution yields the same density as the actual distri- 
bution, n = J2 P f = S P / • From Eq. the deviation 
from the density can be expressed as 



Sn(q,t) = fj (£, lu — 0)5p(q,t), 



(5) 



where n RPA (t, ui = 0) is the polarization in random- 
phase approximation (RPA) at zero frequency with ma- 
terial parameters like electron density and the effective 
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temperature fixed to the value at time t. Note that this 
RPA-polarization enters only the right hand side of QJ. 
The actual polarization which we will derive from is 
more complex. 

The induced potential is given by the Poisson equation 



V ind (q,t) = (p + - 1) 



eo<r 



7 Sn(q,t). (6) 



Now we can solve the kinetic equation for a small 
external perturbation V cxt . We split the distribution into 
the homogeneous part and a small inhomogeneous per- 
turbation, f(p,q,t) = f p (t) + 5f(p,q 7 t). Neglecting the 
time derivative of f p compared to Sf we can solve Eq. 
up to linear order in V cxt 



6f ^ t] = i I dt ' ^ 
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where the perturbation starts at Sf(p, q,to) = 0. 
The density response, x(t,t'), defined by 



6n(q,t)= / dt' x (t,t')V q ext (t') 



(8) 



follows from Eq. J7J by integrating over p with the result 



x (t,t')=Tl(t,t')+ dt 



U(t,t) — + I{t,t) 
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Here Mcrmin's correction is represented by the term 



I(t,t' 



TI RPA (i',0) 
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The polarization, 

p 

describes the response of the system with respect to the 
induced field in contrast to the response function (jHJ) 
which describes the response of the system with respect 
to the external field. 

To link the derived formulas with familiar results we 
note that for time-independent f p , the Fourier transform 
of the polarization with respect to the difference time 

fi(T,w)= fdrre iuT U(T+^, T-~) (12) 



becomes the standard RPA polarization with the fre- 
quency argument shifted by the inverse relaxation time 



TL(T,w) = IT 



RPA 
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Moreover, in the limit t + t' — » oo the solution of Eq. 
© approaches the familiar form \ — n M /(l — ^yiI M ) 
where the Mermin-Das polarization reads 
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Eq. |J[J} describes how the Mermin susceptibility is 
formed after a fast release of free carriers. 

In the further analysis we will follow closely the exper- 
imental way of analyzing the two-time response function 
[3j. The pump pulse is creating charge carriers in the 
conduction band at time to and the probe pulse is sent 
a time to later. The time delay after this probe pulse 
T = t — tp — to is then Fourier transformed into frequency. 
Of course everything starts at to when the pump pulse 
has created the carriers in the conduction band, i.e. the 
maximum delay is T < t — t . The frequency-dependent 
inverse dielectric function associated with the actual time 
t thus reads 
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dTe iujT X (t,t-T). (15) 



This is exactly the one-sided Fourier transform intro- 
duced in Ref. |l5|. It is worth noting that Wigner's form 
of the Fourier transform l|12|) would just result in a factor 
of two in the time evolution of the dielectric function. 

Now we turn our attention to short times when the 
response is built up. In general the integral equation @ 
together with Eqs. ifTUj) and i|TT|) has to be solved. The 
useful limit of long wave lengths q — > offers an appre- 
ciable simplification in that the leading terms of (|11|) and 

JTSJ are U(t, t') « -t)e^ and J(M') w ±e^. 

16] It is practical to transform the integral equation © 
for the response function directly into one for the inverse 
dielectric function, Eq. I|15|l > by differentiating twice with 
respect to time 

z"{t) + (--2iu)z'(t)+\uil(t)-ui 2 -i-} z(t) = -w?(t) 



z{to) = 0; 
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= 0. 
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Here z'(t) — dz(t)/dt and we have denoted the time- 
dependent plasma frequency by ui p {t) = ■ 

In our simplified treatment the density of electrons n(t) 
is assumed to be a known function of the time. Due to 
Rabi oscillations it is a nontrivial function of the laser 
pulse. If the pulse is short compared to the formation 
of the collective mode, all details of the density build-up 
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become unimportant and we can approximate the time 
dependence of the density by a step function 0] . Then 
the plasma frequency is uj p (t) = u> p <d(t — to), where to 
is the time when the pulse reaches its maximum. For 
this approximation the differential equation Ijltil) can be 
solved analytically, yielding 



e(u,t) 



= 1 - 



Y 



dt'e^-^*' sin(rt'), (17) 



where T = yup 1 — 472- The integral l|17|) can be ex- 
pressed in terms of elementary functions, but we find the 
integral form more transparent. 

Formula (|17|l has been derived from Mermin's density- 
conserving approximation. The long-time limit yields the 
Drude formula 



lim ~ = 1 - ~2 7 T-v 



(18) 



This obvious limit is not so easy to achieve within short- 
time expansions. For example the approximate result 
from quantum kinetic theory presented in |l5| gives the 
long-time limit of the form 1 — uj 2 / [ui 2 — (ui + i/r) 2 }. 

In spite of the numerous approximations, Eq. (fTTfl fits 
well the experimental data, as shown in Fig. ^ f° r the 
polar semiconductor GaAs. 

The resonance visible at about 8 THz is due to optical 
phonons. This feature is easily accounted for, see Fig. ^ 
by adding the intrinsic contribution of the crystal lattice 
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eGaAs(w,£) 
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with the longitudinal and transversal optical frequencies 
wlo/2tt = 8.8 THz and wto/2tt = 8.1 THz, the lattice 
damping 7 =0.2 ps" 1 and the non-resonant background 
polarizability of the ion lattice f M = 11.0, see Ref. 0- 

The formula l|17l) results in a too fast build-up of the 
collective mode. This is seen in the second row of Fig. Q] 
at the time t — 2h fs. This is, however, just the time 
duration of the experimental pulse and consequently the 
time of populating the conduction band. Since we have 
approximated this by the instant jump, this discrepancy 
at short times around 25 fs can be expected. A more re- 
alistic smooth populating can be modeled by an arctan- 
function and then the numerical solution of (|16|l improves 
the description as we see in Fig. Only the early de- 
velopment is plotted since the later stages agree with the 
simple estimate <jl"7)) . 

It is obvious that only the gross feature of formation 
of collective modes can be described by the mean-field 
fluctuations. For correlations beyond the mean field we 
should observe deviations. This is the case in the recently 
observed response of InP [3| , where a coupling between 
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FIG. 1: The time evolution of the inverse dielectric function 
in GaAs. The labels from top to bottom denote the time 
t. The pumb pulse was at to = — 40fs and the probe pulse 
has a full- width at half-maximum of 27 fs. Circles are data 
from 0,0, solid lines show the electronic part with intrinsic 
polarization 11911 and dotted lines the electronic part alone 
(E|. The plasma frequency is given by cu p =14.4 THz and 
the relaxation time is r = 85 fs. 



LO phonon and photon modes is reported. In Fig. [3] we 
see that the formation of the collective mode is delayed 
in the experiment during the first 100 fs when compared 
to the mean-field formula H17|) . Such behavior is due 
to scattering processes neglected here. At times 100 fs- 
200 fs the formation of collective modes is described again 
sufficiently well. For InP we use parameters according 
to namely w LO /27r = 10.3 THz, u TO /2ir = 9 THz, 
7 = 0.2 ps _1 , and £00 = 9.6. 

The aim of the present paper was to separate the gross 
feature of the formation of collective modes at transient 
times which are due to simple mean-field fluctuations. 
This has resulted in a simple analytic formula for the 
time dependence of the dielectric function. Subtracting 
this gross feature from the data allows one to extract 
the effects which are due to higher-order correlations and 
which have to be simulated by quantum kinetic theory 
0, S El and response functions with approximations 
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FIG. 2: The first two time evolution steps of Fig. ^comparing 
the instant population of the conduction band (thin line) with 
a smooth transition (solid line). 



beyond the mean field |TjJ. These treatments are nu- 
merically demanding such that analytic expressions for 
the time dependence of some variables are useful for 
controlling the numerics. 

To conclude, we have derived a simple analytic formula 
for the formation of a collective mode. Being able to de- 
scribe universal features of the formation of quasiparti- 
cles, the simplicity of the presented result is extremely 
practical and offers a wide range of applications. We see 
two prominent examples: (1) It could spare a lot of com- 
putational power to simulate ultrashort-time behavior of 
new nano-devices and (2) it can help to understand and 
describe the formation of collective modes during nuclear 
collisions which are not experimentally accessible in the 
early phase of collision. 

We thank A. Huber and R. Leitenstorfer for provid- 
ing us with the experimental data and H. N. Kwong for 
enlightening and clarifying discussions. 
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